Rectification of Fluctuations in an Underdamped Ratchet 
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We investigate analytically the motion of underdamped particles subject to a deterministic peri- 
odic potential and a periodic temperature. Despite the fact that an underamped particle experiences 
the temperature oscillation many times in its escape out of a well and in its motion along the po- 
tential, a net directed current linear in the friction constant is found. If both the potential and the 
temperature modulation are sinusoidal with a phase lag S, this current is proportional to sin 8. 



PACS numbers: 05.40.+j, 82.20.Mj, 87.10.+e 

Recently, research on molecular motors Jl|,||, which 
apparently operate with an efficiency close to kT, has 
stimulated interest in ratchets. Ratchets are commonly 
defined as systems which, in the absence of a net force 
or macroscopic gradient, are able to produce a directed 
current through the rectification of noise j^] . Fundamen- 
tally ratchets are of interest as simple non-equilibrium 
systems. Some of the proposed models of ratchets con- 
sider particles in a periodic potential and subject to a 
non-uniform periodic temperature profile [^-0, or sub- 
ject to a time-dependent external force with zero average 
which, in particular, can be stochastic. Other mod- 
els invoke Brownian motion in fluctuating potentials || 
or models of Brownian particles which can be in sev- 
eral states with external pumping of particles to one of 
the states |l(J. These works treat overdamped particles. 
Models in which inertia plays a role either within classical 
or quantum dynamics have been considered only recently 
0,|2|, and the results available by now are mostly nu- 
merical. 

We present analytical results for the motion of under- 
damped particles in a periodic potential and subject to a 
periodic temperature modulation (see Fig. 1). In the par- 
ticular case shown in Fig. la, both the temperature mod- 
ulation and the potential are sinusoidal with the same 
microscopic period and a relative phase lag S. We show 
that a directed current results which is proportional to 
sin S. No current is observed if the system is symmetric 
5 = mr. This is the underdamped analog of a similar ef- 
fect which occurs for overdamped motion in systems with 
state dependent diffusion . In the overdamped case, 
directed motion appears because a particle diffusing up 
a potential hill can easier surmount a potential barrier 
if the ascent occurs in a heated region |13[. In contrast, 
in the presence of a very small frictional force, a parti- 
cle experiences the microscopic temperature modulation 
many times, both when the particle is in a well and even 
after its escape into a running state. Thus, the current 
must vanish in the zero dissipation limit. 

The motion of a classical particle in an external poten- 
tial and subject to thermal noise is governed by Kramers' 
equation Oj. While the overdamped motion (Smolu- 



chowski limit) is relatively easy to handle, the under- 
damped limit is much more complex. The steady state 
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FIG. 1. One period of the potential U(x), (a) the potential 
given by Eq.(|l7|) (bold solid line), temperature variation T(x) 
(dash-dotted line); (b) an example of sawtooth potential. 



current of particles in a periodic potential due to a con- 
stant external force was treated by Risken and Vollmer 
M|. Here we generalize their approach to the case of 
a spatially non-uniform temperature, thus presenting a 
rare example of a solution of an underdamped problem. 
Kramers' equation Jlif for the joint distribution function 
W(x, v) of coordinate and velocity for a particle with 
mass m in the potential U (x) subject to a frictional force 
with damping constant 7 and a temperature T{x) is 
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= 0. 



(1) 



We assume that both U(x) and T(x) are periodic func- 
tions with the same period L, and the potential U(x) has 
only one maximum per oscillation period. The amplitude 
of the potential is Eq, the zero of energy is chosen at the 
minima of the potential, and the points x = ±L/2 cor- 
respond to the potential maxima (Fig. 1). Furthermore 
we assume that T(x) >C Eq for any x. 

Eq. (Q) must be supplemented by boundary condi- 
tions. For E > Eq the distribution function is periodic 
in x, W(—L/2,v) = W(L/2,v). For energies below Eq 
all particles are reflected from the potential barrier at the 
points x\(E) and x 2 (E) (Fig. 1). Thus, for E < E , the 
boundary condition requires that at the turning points 
the distribution function of incoming particles is the same 
as that of the reflected particles. 

Following Ref. |I7J], we replace the velocity v by the 
energy variable E = mv 2 /2 + U(x). It is then neces- 
sary to discriminate between left- and right-moving par- 
ticles; their distribution functions are denoted as W~ 
and W + , respectively. Introducing for further conve- 
nience the symmetric and antisymmetric combinations, 
W S.A = w+ ±W~,we obtain |J| 



dW A d 

— W 

dw s d 

7— i v{x,E) 



1 + T ^ 



dx 



dE 



1 + T{x) 



dE 



with the boundary conditions 



W s - A (-L/2, E) = W S ' A {L/2, E) , E > Eq 
W A (x 1 {E),E) = W A (x 2 (E),E) = , E<E - 



W s (x,E)j; (2) 
W A (x,E)\, (3) 



(4) 



Here v(x,E) — [(2/m)(E — U(x))} 1 ^ 2 , and the turning 
points x± t 2(E) are the solutions of the equation U(x) = 
E. In addition, both functions W s and W A and their en- 
ergy derivatives are continuous for E = Eq. The function 
W is normalized according to 



dEdx[mv(x,E)]~ 1 W s (x,E) = 1, 



(5) 



and the expression for the (particle) current has the form 

/ = m- 1 J dEW A (x,E). (6) 

Integrating Eq. (^|) over energy, we conclude immedi- 
ately that the current I does not depend on x. Note that 
no assumption about damping has been made so far — 
equations (|2|), (|^) are valid for arbitrary 7. 

Now we turn to the underdamped case. For 7 = 0, 
both the symmetric and the antisymmetric distribution 



functions W s and W A depend on energy only. Then due 
to the boundary conditions W A vanishes for E < Eq, 
and the matching conditions for E = Eq imply that it 
vanishes also for E > Eq . Thus in the absence of friction 
there is no current jl6) . Dissipation is essential for the mi- 
gration of particles between different potential wells. For 
low dissipation the typical time for this migration is much 
longer than the period of the oscillations inside each well, 
and thus the energy varies much slower in time than the 
coordinate x of the particle. Consequently, the distribu- 
tion function depends weakly on x while it varies rapidly 
with energy. This rapid dependence is, however, essential 
only in the narrow layer of energies around E = Eq\ the 
thickness of this layer is proportional to a positive power 
of 7. Generalizing the approach of Ref. JT7J we write 



W S ' A (x, E) = W^ A {x, E) + w^ix, E) 



(7) 



where W S ' A are slow functions of x and E, while 
w S ' A (x, E) are slow functions of x that vary rapidly with 
energy. The functions w S ' A are different from zero only 
in the narrow layer of energies around Eq ; their role is to 
ensure the continuity of W S ' A for E = E . For E > E , 
the boundary conditions read now 



W ' (-L/2, E) = W b ' A (L/2, E) 



w 



S,A 



(-L/2,E) 



,.S,A 



(L/2,E), 



(8) 



and for E < Eq 

W A {x ia (E), E) = w A ( Xl a{E), E) = 0. (9) 

First we consider the functions W S ' A (x, E). The pe- 
riodic functions v(x,E), t(x,E) = v(x, E)T(x) and the 
distributions W S ' (x, E) can be expanded in Fourier se- 
ries (in the range E < U{x) we define them to be zero), 



00 ^ 00 

v(x, /•.') c, • ^ v nc cos — - ^ 



L ^ 
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W s ' A (x, E) = W^ A {E) + £ W*> A (E) cos 



2irnx 



sin 



2ttux 



It is seen from Eqs. (^), (^) that the coefficients W A C s are 
proportional to 7. Then for E < Eq due to the bound- 
ary condition the whole function W A is proportional to 
7, and consequently the x-dependent part of W s is pro- 
portional to 7 2 . In the following we restrict our consider- 
ation to the terms proportional to 7, thus neglecting the 
coefficients W^ c s - Taking the constant part of Eq. (j^), 
we find for E < Eq 



2 



W S (E) ee NP(E) = iVcxp { - f E j^dE' 



(10) 



with a normalizing constant N. With the abbreviation 
(f>(x,E) = 1 - (v (E)/t (E))T(x) the solution for W A 
satisfying the boundary condition for E < Eq is 



W A (x,E)= 7 JLlw s (E) 



xi(E) 



v(x',E)(j){x',E)dx l 



We note that for potentials U(x) (such as that of Fig. la) 
with a saddle at Eq, this expression diverges as E — » Eq 
and x — > ±I//2. This divergence is similar to that of 
the oscillation period for E Eq. It is not detrimental, 
since the Jacobian from (x,v) to the (x,E) space has a 
compensating singularity. It is advantageous to treat this 
divergence by considering a sawtooth potential, like that 
of Fig. lb, for which all expression converge. We will pro- 
ceed in this way, and our final results will be shown to 
be convergent in general case. Thus, we expect that our 
consideration is also valid for potentials with a saddle. 

For E > E Q the function W still does not depend on 
x (up to terms of order 7), and is given by Eq. ([To]). For 
W A [x,E) we find 



W A (x,E)= 1 -^^W s (E) 
+ B(E), 



L/2 



v(x',E)(f)(x',E)dx' 



(11) 



which differs from the solution for E < Eq by the 
position-independent function B(E). The latter appears 



due to the different boundary conditions, and is deter- 
mined by Eq. (||). Taking the constant part of this 
equation and fixing the only integration constant so that 
B(Eq) — 0, we obtain 



B(E) = P(E) 



F(E 



J-dE 1 + G{Eq) 



>e P(E> 

Here we have introduced the functions 



-G(E). (12) 



1 OO 

x ' n— 1 %— c,s 



8W A 

v ni (E)W A +t ni (E)-^f 



G(E) 



Llky^ J ^ /2 xv (x,E)0(x,E)dx 



The functions W A C s (which are just the coefficients in 
the Fourier series for the solution Eq. (|ll|)) are easily 
obtained from Eq. (0) and read 



w A =±lL± 

nc ' s 2m dE 



d 

v ns , c (E) + t nSiC (E) — 



W S (E) 



With the help of the antisymmetric product, 

r L/2 



rL/2 

{a x b} = (2L 2 )^ 1 / a(x)b{x')sign{x — x')dxdx ', 

J -L/2 



the expression for B(E) after some algebra can be 
brought into the form 



B(E) = -jLNP(E) 



E 



dE' 
toW) 
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(13) 



Now we cure the discontinuity of 8W A /dE for E = Eq by 
taking into account the functions w S ' (x, E). Since they 
vary rapidly with energy, the derivatives dw/dE can be 
eliminated in favor of d 2 w/dE 2 . Replacing E by Eq in 
the functions which are smooth in the energy, neglecting 
the difference between X\{E) and —L/2, and introducing 
the dimensionless variable it, 



u(x) 



2tt 



I-l/2 T ( x ) v ( x i E )dx 



(which is an analog of the action for the uniform temper- 
ature case), — 7r < u < 7r, we find 



dw S ' A 
du 



2„„A,S 



d w 



Uq- 



dE 2 



(14) 



Here uq = r yLto(Eo)/(27v), and the boundary conditions 
read w a (±tt, E) — for E < Eq, and w s > a (—tt,E) = 



w 



S,A 



(tt,E) for E > E . 



Eq. (|14| ) has been investigated by Risken and Vollmcr 
p5| . Making use of their solution, writing down the 
matching conditions for W A at E — Eq, and taking into 
account that the function B(E) in Eq. ( |ll| ) is propor- 
tional to 7, we find that the functions w • are propor- 
tional to 7 3 / 2 . Thus, these functions are not essential for 
our consideration. For the current we obtain then 



I = m - L I B(E)dE. 

IE 



(15) 



Eqs. @ and © 

together with the definition (g|) of the 
normalizing constant N yield a closed expression for the 
current induced by an arbitrary temperature distribution 
T(x). This current is proportional to 7 and vanishes in 
the limit of constant temperature T(x) — Tq, since then 
t{x,E) = Tqv(x,E). 

We continue with two simplifying assumptions. First, 
we consider the experimentally relevant regime of small 
temperature gradients. Specifically, following Rcf. PJ, we 



3 



assume 



T(x) =T + T lC os(2irx/L~6), T x < T . 

Furthermore, for simplicity, we choose the potential U (x) 
to be symmetric, U{x) — U{—x). To linear order in Zi, 
the expression for the current simplifies and is given by 



/ = 



jLNTi 



mT L 2 v (E ) 
v{x' , Ef))dx' 



sin 5 



f L/2 ■ I 




/ dx sin 


^2ttx^ 


l-L/2 





dE 



(16) 



where we used the condition Tq <C Eg. As is expected, 
the current is proportional to sin (5, and vanishes when 
both the potential and the temperature are symmetric. 
We also note that Eq. ([f6j) is finite for a potential U(x) 
of a general form (even though dv(x, E )/dx diverges for 
x — > ±i/2), and thus, we expect that it is generally valid. 
As an example, we consider the simple form of U(x), 



U(x) = (E Q /2) [1 - cos(27nc/L)] , 
shown in Fig. la. We obtain 



(17) 



2mTo 



exp 



-— I Mil-). 



(18) 



It is instructive to compare this expression with the one 
found in the overdamped regime f^fy (index ov). For the 
same parameters we find 



^ElT x 



exp 



sin i5. 



(19) 



First, we note that the sign of the current is the same 
in the overdamped and in the underdamped regime. It 
corresponds to the notion jl3],Q that it is easier for the 
particle to climb a "hot" slope. Furthermore, we note 
that the two expressions fll8| ) and (|ll]) match for a cer- 
tain critical damping, 

7c ~ (E /L)(m/T )^ 2 . 



This critical damping can be understood by considering 
the overdamped side. Indeed, to leading order, the over- 
damped approximation neglects the terms which 
do not contain 7 in Eq. ([!]). This is justified for 
7 3> v~ 2 dU/dx, where v ~ (T/m) 1 / 2 is the typical ther- 
mal velocity. Substituting Eq/L for the characteristic 
value of the potential derivative, we readily obtain Eq. 
©. Thus, we conclude that our expression for the un- 
derdamped regime (^) is valid up to 7 ~ j c . 

Now we compare the result ( |l8| ) with the current 
I ~ (F/jL) exp(-E /T ), obtained in Ref. (l5) for the 
underdamped Brownian motion due to an external force 
F. We see that a non-uniform temperature has the 
same effect as an applied effective uniform driving force 



F e ft ~ j 2 LTi sin S/(mTo). Thus, the motion generated 
by a non-uniform temperature profile is less effective than 
the forced motion due to a uniform field. This is in con- 
trast to the overdamped case Q , where the effective force 
F e ff ~ E0T1 sin 5 /(T L) does not depend on 7. 

In conclusion, we investigated how a non-uniform tem- 
perature can produce a steady, directed current of parti- 
cles in a periodic potential in the underdamped regime. 
A current is generated even in the case of a symmetric 
potential provided the temperature profile is asymmetric 
(out of phase). A temperature oscillation is equivalent to 
an external driving force with a magnitude proportional 
to the square of the friction constant. 

The work was supported by the Swiss National Science 
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